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Abstract. A physical interpretation of the rope simulated by the 
SONO algorithm is presented. Properties of the tight polygonal 
knots delivered by the algorithm are analyzed. An algorithm for 
bounding the ropelength of a smooth inscribed knot is shown. Two 
ways of calculating the ropelength of tight polygonal knots are 
compared. An analytical calculation performed for a model knot 
shows that an appropriately weighted average should provide a 
good estimation of the minimum ropelength for relatively small 
numbers of edges. 



1. Introduction 

Knots tied on material objects, e.g. a rope of a finite thickness, are 
called physical knots JB]- It is well known from everyday life that 
open knots can be tightened by pulling apart the loose ends of the 
rope. Tightening a closed knot, i.e. a knot tied on a rope whose 
ends are spliced, is more difficult since one has to engage processes 
in which the rope will reduce its length while keeping its diameter 
intact. Such a process, difficult to create in the laboratory, can be easily 
simulated numerically and we find in the literature various examples 
of such simulations [IB H d d 113 EH El ■ From the intuitive point 
of view, it is rather obvious that in tightening a knot tied on a rope 
of a finite diameter, one must arrive at a conformation at which the 
tightening process stops. This conformation is called tight. Any small 
change in the shape of a tight conformation needs an elongation of the 
rope. Obviously, tight conformations of knots tied on different types 
of rope, with different physical properties, will be different. Thus, to 
make the notion of the tight knot unambiguous, one must specify the 
physical model of the rope on which the knot is tied. The simplest of 
such models is the perfect rope. Perfect rope is perfectly floppy - needs 
no force to be bent, but, at the same time, it is perfectly hard - needs 
an infinite force to be squeezed, its diameter always remains intact. It 
is also perfectly slippery. Its formal definition will be given in the next 
section. 

l 
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Tight conformations can be seen as minimizers of the ropelength 
function [11 . The minima can be only local. While it has been 
shown that global minima exist [3J E] , with the exception of the triv- 
ial knot and some simple links, the minimizing conformations are not 
known. Knots in the conformations at which the ropelength reaches 
its global minimum are called ideal. As long as the ideal conformations 
are not known analytically, we cannot be certain that the conformation 
at which one arrives is ideal. 

In all simulations of the processes in which knots are tightened, the 
knots are represented by a finite sequence of points in M. 3 space. When 
connected, the points can be seen as vertices of a self-avoiding polygon. 
Most often the simulation programs either start from a non-equilateral 
conformation and aim at making it equilateral (e.g. SONO) or start 
from an equilateral conformation and try to keep it such (crankshaft 
rotations). Let P be such an equilateral polygon. The essential prob- 
lem we face is to construct, using the coordinates of the vertices of 
P, a smooth knot Kp which can be seen as the axis of the knot tied 
on the perfect rope. In constructing Kp, we make some assumptions 
concerning the shape of the pieces of which it is built. In what follows 
we assume that the knot Kp is built from smoothly connected arcs in- 
scribed into P, which makes the knot C 1 ' 1 smooth. As a result of this 
assumption and provided the radius of the rope is properly chosen, we 
show that the knot Kp can indeed be tied with perfect rope leaving it 
self-avoiding. On the other hand, we also know that the inscribed knot 
is not ideal, thus, its ropelength can serve only as an upper bound for 
the minimum ropelength of the unknown ideal conformation. In find- 
ing tight polygonal knots P with an increasing number of vertices, we 
see that the ropelength of the inscribed smooth knots Kp converges to 
a value, which we believe is the ropelength of the ideal conformation, 
although we will not prove that here. This convergence is slow however, 
and one might want an approximation using relatively few edges. This 
is especially important for complicated knots where the number of edges 
needed to adequately approximate the smooth ropelength would not be 
practical. With this aim, a second, weighted-average, approach is ex- 
plored. This technique provides surprisingly accurate approximations 
with relatively few edges. Together, we provide both a theoretically- 
sound approach with slow convergence and a practical approach with 
fast convergence. Depending on the application, we expect that both 
techniques will be useful. 
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2. Computing upper bounds for smooth knots inscribed in 

tight polygons 

2.1. Ropelength of knots tied on the perfect rope. Smooth rope- 
length models perfect rope as a non self-intersecting tube with a C 1 ' 1 
knot as its core. For completeness, we include these definitions and 
some properties of this knot energy. 

Definition 1. For a C 1 knot K and x G K, let D r (x) be the disk of 
radius r centered at x lying in the plane normal to the tangent vector 
at x. Let 

R(K) = sup{r > : D r (x) n D r (y) = for all x ^ y e K} . 

The quantity R(K) is called the thickness radius of K. Define the 
ropelength of K to be 

Rope(K) = Length(K)/R(K) , 

where Length(K) is the arclength of K. 

The thickness radius R{K) is the radius of a thickest tube that can 
be placed about the knotted core K without self-intersection. Note 
that in some of the literature, the thickness refers to the diameter of 
this thickest tube, in which case the ropelength is half of what we use 
here. The problem we face is how, knowing the shape of K, to find 
R{K). 

Suppose we are given a concrete knot K and we are trying to repro- 
duce it with an impenetrable tube. Thus, K serves here as the core of 
the knotted tube. We have to pay attention to two types of interactions 
between the tube's normal disks that restrict the radius R of the tube. 

(A) The interactions between the discs centered on pairs of points 
located infinitesimally close to the diagonal of K x K - these 
normal disks will start overlapping when the curvature radius 
is smaller than R. 

(B) The interactions between the discs centered on points located 
away from the diagonal of K x K - the disks will start overlap 
when the Euclidean distance between the points is smaller than 
2R. 

Figure ^ illustrates the problem. These intuitive observations are 
captured by the quantities below and the subsequent lemma. 

Definition 2. For a C 2 knot K with unit tangent map T, let MinRad(K) 
be the minimum of the radius of curvature measured at each of the 
points on the knot. The doubly- critical self-distance is the minimum 
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FIGURE 1. Two different conformations of the trefoil 
knot are inflated. Both conformations have the same 
arclength. The inflation of the conformation shown on 
the left stops first because of its sharp bends at which 
interactions of type A intervene. Further inflation would 
lead to the development of singularities in the tube sur- 
face. See the framed picture. Inflation of the conforma- 
tion shown on the right stops at the point above which 
interactions of type B intervene. 

distance between pairs of points on the knot whose chord is perpendic- 
ular to the tangent vectors at the both of the points. In other words, 
let 

DC{K) = {(x,y)eK xK : T(x) ±xy± T(y),x ^ y} 

where xy is the chord connecting x and y. Define the doubly- critical 
self- distance by 

DCSD(K) = mm{\x-y\ : (x, y) G DC(K)}, 
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where | • | is the standard IR 3 norm. 

There is a fundamental relationship between R(K), MinRad(K), 
and DCSD(K). 

Lemma 3. Suppose K is a C 2 knot. Then 

R(K) = min {MinRad(K), DCS ^ K ) j . 



Proof. See HU. □ 

In jSJUj, it is shown that ropelength minima exist as (at worst) C 1 ' 
curves. 

After reworking the definition of MinRad, one can extend Lemma |3] 
to include C 1 (and thus C ' l curves). The following is taken from j^j. 
For a C° function / : R — > M 3 , define the dilation of / by 

dil(f) =sup( l/( ^ : s,ieR,^t 

I I s ~~ *l 

Note that if K is a C 2 knot parameterized by arclength with unit tan- 
gent map T, then MinRad(K) = l/dil(T). Thus, the dilation gives a 
generalization for MinRad to knots that are C 1 . Since MinRad and 
l/dil(T) are equal for C 2 knots, we use MinRad to denote l/dil(T) 
for all C 1 knots. If a knot is C 1 but not C 1,1 , then dil(T) = oo, in 
which case MinRad is assumed to be 0. In this paper, we are mainly 
interested in C 1,1 knots, in which case dil(T) is finite and MinRad is 
positive. Alternate, but equivalent, approaches for defining the rope- 
length of C 1,1 knots are explored in 012] • 

Lemma 4. Suppose K is a C 1 knot. Then 

R{K) = min {MinRad{K), DCS ^ \ K > \ 



Furthermore, if K is C 1 ' , then R{K) > 0. 

Proof See |3] or [§]. □ 

2.2. Ropelength of polygonal knots. As mentioned in the intro- 
duction, the problem we face in analyzing knots provided by numerical 
simulations is that they are not given in a smooth continuous form. For 
instance, programs based on the SONO algorithm deliver equilateral 
polygonal knots. Let P be such a knot. It is represented by a sequence 
of n points whose positions are indicated by vectors v±, V2, ■ ■ ■ , v n . We 
shall refer to them as the vertices or beads of the polygonal knot P 
that one obtains by joining consecutive vertices. We define 
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Si 


= V i+ i - Vi 






Si 


= \Si\ 






r i,j 


= Vj - Vi 






T i,j 


= \ r i,j\ ■ 



The vectors will be called edges. When referring to the vertices as 
points on a polygonal knot, we shall omit the vector notation. 
Remark: Because the knots that we are considering are closed, when- 
ever the vertex index k happens to be larger than n, we assume that 
Vk = Vk-n- Similarly, when the index happens to be smaller than 1, we 
assume Vf- = Vk+ n - 

In what follows, we will use the notion of the index distance ID 
between the vertices: 

% — j\ , when \i — j\ < n/2 



TT^t \ f \ l — J ; when \i - 
v ■" I n— \i — j\, otherwise 



Note that alternate approaches have been appeared in 12U E3 IHB 
Uni 13 El HZl Uni HI In the SONO simulations, the vertices of a knot 
interact with each other as if they were spheres, the properties of which 
will be later. Taking into account the representation, we shall extend 
the notion of the thickness radius in terms of the virtual spheres. We do 
this by defining notions of MinRad and DCSD for polygons similar to 
the characterization theorem (Lemma HJ for smooth thickness radius. 

Definition 5. Suppose P is an equilateral polygonal knot with edge 
lengths dl. Let 

Rad{Vi) = 2tan(U) 
where 0, is the turning angle of the tangent vectors at Vi. Note that 
Rad(vi) is the radius of an arc of a circle a that can be inscribed in the 
corner at Vi so that a intersects both adjacent edges at the midpoint 
and the tangent at the midpoints coincide with the tangents on a at 
the midpoints. Let 

MinRad(P) = mm{Rad(v {)} = — — , 

i 2tan(0 maa; /2) 

where 9 max is the maximum turning angle. See Figure El 

Definition 6. Let P be a polygonal knot. We define the sphere thick- 
ness radius to be 

R 8 (P) = min{V 'MinRad(P) 2 + (dl/2) 2 , SR(P)} 
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Figure 2. Vertices and edges of a piece of an equilateral 
polygonal knot P. Inscribed arcs of radii Rad(vi) are also 
shown. To simplify the drawing we assumed that all of 
the vertices are located in the same plane. 



where SR(P), the sphere radius, is half the minimum distance between 
vertices whose index distance is at least \irMinRad(P)/dl — 1]. 

The construction of the sphere thickness radius needs a few words 
of explanation. Suppose P is an equilateral knot. Imagine that about 
each of the vertices, there is a sphere of radius R. When R is larger than 
dl/2, the spheres about nearby vertices intersect. The cases in which 
dl/2 < R < dl/y2 need special attention, but since they never occur 
in the simulations we eventually analyze, we exclude them assuming 
implicitly in what follows that R > dl/y/2. The surface of the union of 
all spheres is bumpy, so we refer to this as the corrugated tube or bead 
rope about the knot. The problem we face is to inflate the spheres as 
much as possible without violating the condition of the self-avoidedness 
of the corrugated tube formed by their union. 

First, let us explain what we mean by self- avoiding in the case of 
the corrugated tube. Consider three consecutive vertices v i: and 
Vi+i and the spheres 5i_i, Si, and Si+i, around them all of radius R. 
The spheres Si-i and Si intersect at a circle Cj_i whose radius equals 
Rb = a/-R 2 — {dl/2) 2 . The spheres Si and S i+ \ intersect at a circle Cj 
of the same radius. Let -Dj_i and Di be the disks filling the circles. 
The disks are inclined at an angle equal to the turning angle di. We 
assume that the disks are hard and hence cannot overlap; thus, at most 
they may become tangent at a single point. The latter happens when 
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Figure 3. Violation of the self-avoidedness of the bead 
tube, (a) The tube bends too quickly, (b) Two index- 
distant portions of the tube are too close to each other. 



the sphere radius Ri = 2s m(6> /2) • ^ or a gi ven polygonal knot P, the 
smallest of such radii is y/ MinRad(P) 2 + (dl/2) 2 . When R is chosen 
larger than ^ MinRad(P) 2 + (dl/2) 2 , we see the behavior shown in 
Figure Et- 

When R < yj ' MinRad(P) 2 + (dl/2) 2 , all portions of the corrugated 
tube are well-defined locally and one must check if index-distant ver- 
tices are too Euclidean-close. Let us explain, what we mean by stating 
that the portions of the corrugated tube are well-defined. Each por- 
tion of the corrugated tube is a sphere with two pieces cut off. The 
cutting planes are perpendicular to the polygon edges and are located 
at the midpoints. Consecutive portions are thus connected with disks 
of radius R^. As stated above, the disks are assumed to be hard and 
they are not allowed to overlap. A portion of the corrugated tube 
which fulfills the condition is said to be well-defined. If the Euclidean 
distance between two vertices is smaller than 2R, the portions of the 
corrugated tube would overlap. Thus, we check all the Euclidean dis- 
tances between sufficiently index-distant vertices and find the smallest 
one. This will be the value of 2SR(P). This is a geometrical interpreta- 
tion of the self-avoidedness conditions found in the definition of R S (P). 
When R is chosen larger then SR(P), we see the behavior shown in 
Figure Eb- 

We will be mainly interested in optimizing R s , but to do so we need 
a scale-invariant version of R s . There are various ways to do this and 
we explore two normalizations in the next section. The corrugation 
will play a part in how the knot packs efficiently in these optimizing 
conformations. In the short term, our goal is to show that for a suffi- 
ciently "thick" polygon P, there exists a smooth knot inscribed in P 
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whose thickness radius is close to R S (P). We first present the inscrib- 
ing algorithm. A similar formulation for inscribing polygons appears 



Proposition 7. For a given n-edge equilateral polygonal knot P with 
edge length dl, a C 1 ' 1 curve Kp can be inscribed in P in such a way that 
MinRad(Kp) = MinRad(P) . Furthermore, there exists a bijection 
between Kp and P so that for each pair x' G P and x G Kp, we have 



Proof. An arc an of a circle of radius Rad(vi) can be inscribed at t>, such 
that cii is tangent to and s*j and intersects the adjacent edges at 
the midpoints. Let K P be the result of inscribing arcs of radius Rad{vi) 
at each vertex Vi and removing the bypassed corners. Since there is no 
overlapping of adjacent inscribed arcs, K P is well-defined as a (possibly 
self-intersecting) closed curve. The curve Kp has a piece-wise constant 
radius of curvature, and MinRad(Kp) = MinRad(P) . The knot Kp 
is C 1 and piecewise C 2 , and thus, lies in the category of C 1,1 curves. 

For each x on the inscribed curve Kp, we define a unique point x' on 
P. If x is a midpoint on P, let x' = x. Otherwise, x lies on an arc, say 
«i, whose center is Cj. Let x' be the intersection of the ray CiX with 
ej_i U e, (see Figure EJ). Simple trigonometric calculations show that 



For each x G Kp, we say the corresponding vertex is the vertex which 
generates the arc on which x lies. The midpoint of each edge has two 
corresponding vertices. It will be convenient to think of the midpoint 
as corresponding to each vertex in different situations and this will not 
pose any problems in this work. Note that the smaller 8 max is, the 
closer the inscribed curve Kp will be to the polygon P. 

The corrugated tube about the polygon is a union of many spheres, 
some of which are allowed to intersect. The intersection between the 
surface of two adjacent spheres consists of a circle whose radius is 
a/ R s (P) 2 — dl 2 /A and whose center is the midpoint of the edge. This 
value will be important, so we define R C (P) = ^ R S (P) 2 - dl 2 /A. We 
will show now that the smooth tube of radius R C {P) about Kp is non- 
self intersecting, i.e. that the thickness radius R(Kp) of Kp is at least 
R C (P). See Figure El This is a bit delicate and will take some effort to 
prove. We will need the following technical lemma. 



in U7|. 



x-x'\< R S (P) (sec (^) - 1) . 




□ 
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Figure 4. The definition of the points x and x' used in 
Proposition [7| 

Lemma 8. Let P be an equilateral polygonal knot with edge lengths 
dl and R S (P) > dl/\/2 and Kp be the smooth knot inscribed in P 
via the algorithm of Proposition [7| If x,y G K P and arc(x, y) > 
7iMinRad{K P ) = ixMinRad(P), then ID(v x ,v y ) > \nMinRad{P) / dl- 
1], where v x and v y are the vertices corresponding to x and y. 

Proof. Suppose arc(x,y) > nMinRad(P). Now arc(x,y) is the mini- 
mum arclength on Kp measured over the two paths connecting x and 
y. By the construction of the inscribed smooth knot, we know that 
arc(x, y) is no larger than the arclength along P between the corre- 
sponding points x' and y'. In the worst case, both x and y are mid- 
points and their corresponding vertices v x and v y lie in the shorter arc 
between x and y. In such a case, we have 

arc(x, y) < arcp(v x , v y ) + 2(dl/2) , 

where arcp denotes the minimum arclength when walking along the 
polygon P. Now arcp(v x , v y ) = ndl where n is the number of edges be- 
tween v x and v y . Thus, ndl > 7rMinRad(P)—dl or n > nMinRad(P)/dl— 
1. ' □ 

In other words, if arc(x,y) > TcMinRad(P), then we know that the 
corresponding vertices v x and v y will be a part of the set over which 
SR(P) is computed. This will be important in ensuring that SR(P) 
bounds DCSD{K P ). 
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Figure 5. The picture on the left presents the polygonal 
trefoil knot P tightened by SONO. The distance between 
its vertices is guarded by the virtual spheres of radius 
SR(P). Since in the case shown in the figure SR(P) < 
^/MinRad(P) 2 + (dl/2) 2 , the value of SR(P) serves as 
the sphere thickness radius R S (P). Consecutive spheres 
meet at circles of radius R C (P) = \J R S (P) 2 — dl 2 /A. The 
smooth knot Kp shown also in the figure consists of arcs 
inscribed into K. The picture on the right shows the 
inscribed knot Kp inflated to radius R C (P). It is well 
visible, and we prove it in the text, that the smooth 
tube knot, hidden as a whole inside the corrugated tube 
made of the virtual spheres, remains self-avoiding. 



We now have the tools to bound the thickness radius of the inscribed 
knot Kp. 

Theorem 9. Let P be an equilateral polygonal knot with edge lengths dl 
and R S {P) > dlj\[2 and Kp be the smooth knot inscribed in P via the 

algorithm of Proposition^ Then R(K P ) > R C (P) = R S (P) 2 - *f . 

Proof. By the construction of Kp, we know that MinRad(Kp) = 
MinRad(P) . We split the proof into two cases: when DCSD(Kp)/2 > 
MinRad(Kp) and when DCSD(K P )/2 < MinRad(K P ). 

In the former case, we have R(K P ) = MinRad(K P ) = MinRad(P) > 
R C (P). 
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In the latter case, we assume R(Kp) = DCSD(Kp)/2. By [T7j . 
we know that when DCSD(Kp) is realized at a pair of points, say a 
and b, then arc(a,b) > nMinRad(P). Therefore, the corresponding 
vertices v a and have ID(v a ,Vb) > \irMinRad(P)/dl — 1] by Lemma 
IH1 Let D a and D b be the disks of radius R C (P) normal to K P at a 
and b respectively and let B a and B^ be the spheres of radius R S (P) 
centered at v a and v b respectively. Now a and b lie inside B a and 
respectively, and by the definition of R S (P), we know that the interiors 
of B a and B^ do not intersect. Furthermore, the normal disks D a and 
Df, are contained within B a and Bj, respectively. Therefore, D a and 
can only intersect, in the worst case, on the boundary. This implies 
that \a-b\> 2R C (P) or DCSD{K P )/2 > R C (P). □ 

The previous result is not sharp in general. However, for the sake of 
this paper, it does provide a computable lower bound for the thickness 
radius of Kp. Furthermore, the length of Kp, which is smaller than the 
length of P, can be computed explicitly. We call Length(Kp) / R C (P) 
the inscribed ropelength and denote it L c . The value of L c is an upper 
bound for the ropelength of the inscribed Kp, and as dl — > (note that 
changing dl affects n, P, and Kp), our bound on R(K P ) approaches 
R S (P). Since R S (P) > dl/\/2, the inscribed arcs lie within the cor- 
rugated tube about P, which guarantees that Kp has the same knot 
type as P. Since Kp is a C 1,1 smooth knot with the same knot type as 
P, we have a bound on the ropelength of one smooth knot within the 
knot type. Thus, we have an upper bound for the minimum smooth 
ropelength for the given knot type. 

3. Interpretation of simulations performed with the 

sono algorithm 

Let us now describe the SONO algorithm using the notions defined 
above. First notice that in discussing the ropelength of smooth knots, 
we are considering knots K of a given length Length(K) looking for 
their thickness radius R(K). The value R(K) is the maximum ra- 
dius for which K could be inflated without violating the conditions 
that the surface of the tube must remain self-avoiding. The value 
Rope(K) = Length(K) / R(K) is the ropelength of K. The applica- 
tion of the virtual inflation suggests numerical calculations in which 
the inflation process forces the knot to change its conformation. Such 
changes may maximize R{K) and thus minimize Rope(K). Working 
with polygonal knots, it proves more convenient to consider a different 
scheme in the simulations: keep the sphere radius SR fixed and shorten 
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the edge length dl so that the knot eventually arrives at a conformation 
for which R S (P) = SR. This is the approach we discuss below. 

3.1. Basic procedures of SONO. There are two essential proce- 
dures on which SONO is based The basic goal of SONO is to re- 
duce Length(P) subject to the constraint that R S (P) > SR, for some 
fixed rope radius SR. In practice, we choose SR = 1. The reduction 
of the polygon length is achieved by reducing its edge length dl. We 
aim at simulating equilateral knots; thus one of these procedures called 
EqualizeEdges (EE), checks the length of the edges, and, if they differ 
from the desired value dl, introduces necessary corrections. Suppose 
the distance Sj between vertex i and % + 1 is different from the desired 
dl value. Then the vertices are shifted toward new positions so that 
their distance is closer to dl: 

v'j = Vi- c x {dl - Si 
v'i+i = v i+ i + cx(dl - Si) Si/ Si 

c\ G (0, 1/2]. Usually we work with c\ = 1/2. 

As a result of multiple applications of the procedure, the dispersion 
of the edge lengths in the final conformations becomes very small, and 
thus, we consider the knots we get from our simulations as equilateral. 

The second procedure called RemoveOverlaps (RO), checks the dis- 
tances between the vertices. When the RO procedure finds that the 
minimum distance 2SR is violated, the vertices V{ and Vj are shifted 
away from each other to a distance equal to 2SR or, what proves to 
speed up the initial stage of the tightening process, exceeding 2SR by 
e: 

v'i = Vi- c 2 {2SR - r itj + e) r^/r^ 
v 'j = Vj + c 2 {2SR - r id + e) Tij/rij 

where c\ G (0,1/2] and is usually set equal to 1/2. The value of e 
changes during the tightening process. Initially it is of the order 10~ 2 
and reaches the level of 10~ 7 at its end. 

The procedures RO and EE are in some circumstances contradictory, 
but their multiple application leads to simultaneous reduction of both 
the overlaps and the dispersion of the edge lengths. Both parameters 
are constantly monitored. 

When the edge length dl is small, the tightened polygonal knots 
develop short regions where Rad(vi) tends to become small. An addi- 
tional procedure ControlCurvature (CC) monitors this, never allowing 
the external angle between the consecutive edges to be larger than 
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2arcsiii2|^- This, in particular in the final stages of the tightening, 
makes sure that the tight conformation of P will have R S {P) = SR. 

3.2. Physical sense of the SONO algorithm and practical de- 
tails of simulations. Let us now discuss the physical sense of the 
simulation based on the SONO algorithm. The simulated knot P can 
be seen as tied on a closed necklace of beads. The necklace is unusual, 
because its beads are unusual: if their index distance ID(vi,Vj) < 
\ 7TMm ^ ld ( p ^ ~\ , they are allowed to overlap, otherwise they behave as 
hard spheres of radius SR and repel infinitely hard. The RO procedure 
simulates this interaction. Because of the action of the EE procedure, 
the centers of the consecutive beads can be seen as connected with 
nonextensible rods of controlled length dl. The rods are connected to 
each other by elasticity free but angle-limiting joints. Each of the beads 
has the shape of a sphere with two parts of it cut off. The cutting planes 
are located in the middle of the edges which lead to it. Consecutive 
beads are thus connected via disks of radius R c = ^SR 2 -dl 2 /A. The 
necklace can be bent in any direction but only to some extent: the CC 
procedure limits the angle between the consecutive rods to 2 arcsin ^f^. 
As a result, MinRad(P) > R c , which prevents the disks separating 
consecutive beads to overlap. At most they are allowed to become 
tangent. Because of its construction, the surface of the bead rope is 
not smooth. It is corrugated and the corrugation is more pronounced, 
the larger the distance between consecutive beads. This distance is, of 
course, the edge length dl. 

The tightening process runs as follows. Suppose we start with a 
loose conformation. The EE and RO procedures make the lengths of 
all edges equal and remove overlaps between the beads. After multiple 
applications of the procedures, the dispersion of the edge lengths and 
the sum of all overlaps fall below a fixed level. Then SONO reduces 
the required edge length dl. Now, as the EE procedure tries to adjust 
the length of the edges to the new value by shortening the distance 
between consecutive beads, new overlaps may appear. Subsequently 
they are removed by the RO procedure. And so on. 

In the final tight conformation, the dispersion of edge lengths and the 
sum of all overlaps are smaller than 10 -6 which allows us to consider the 
closest beads as just touching. Keeping in mind that the CC procedure 
does not allow a/ MinRad(P) + dl 2 /4 to be smaller than SR, we may 
consider the final tight knot as a polygonal, equilateral knot P whose 
R S {P) = SR. 

The shrinking rope forces the knot to change its conformation. At 
the end, we arrive at a conformation for which further shortening is no 
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FIGURE 6. Tight conformations of the trefoil knot tied 
on a rope consisting of 15, 30, 60, 120, 240 and 480 beads 

longer possible because it creates non-removable overlaps. Obviously, 
if the number of beads is small, the final value of dl is large and the 
bead rope is strongly corrugated. This may create problems, since the 
knot is more likely to become stuck in a local minimum. The reason 
for this is as follows: one part of the rope winding around another part 
may get into a groove between consecutive beads even if shifting it to 
another groove would allow further shortening of the rope. The shifting 
will not be achieved by SONO, because it would need a (rather large) 
temporary increase of the ropelength. Working with a small number of 
beads, we set e to a higher value to minimize these effects. In general, 
the tightening process should be seen as a physical experiment in which 
the experimenter watches and adjusts the parameters. 

3.3. The problem of finding the right ropelength, an experi- 
mental approach. In Figure El we have the tightest conformations of 
the trefoil knot found by SONO working with rope consisting of n =15, 
30, 60, 120, 240, and 480 vertices. 

It seems plausible that the best estimation of the minimal length of 
the rope needed to tie a particular knot will come from analyzing the 
polygonal conformation obtained for highest n. As seen in the figure, 
at n — 480 the surface of the bead rope is visually smooth - the corru- 
gation of its surface, so visible at n = 15, 30, 60, becomes undetectable 
by the naked eye. The problem we face is that in tightening much 
larger knots, for instance the (2, 99) torus knot, reaching this level of 
the surface smoothness would require a very large number of beads, 
which is both awkward and time-consuming. For practical reasons, it 
is reasonable to limit the number of beads as much as possible. We 
present a picture showing the details of the SONO tightened trefoil for 
such a small n in Figure 



16 



JUSTYNA BARANSKA, PIOTR PIERANSKI AND ERIC J. RAWDON 




Figure 7. Details of a trefoil knot tightened by SONO. 
The knot is tied on a piece of corrugated rope consisting 
of n = 15 beads. The centers of the beads make the ver- 
tices of the polygonal knot. The vertices are connected 
with straight cylindrical edges. 

Calculations run faster when n is small, but what about the accuracy 
of the ropelength we can extract from the analysis of the final, tight 
conformation? 

4. Ropelength of SONO knots 

We are mainly interested in understanding the optimal conforma- 
tions of R s . However, R s varies with scale, so one could always increase 
R s simply by scaling the polygon. In the second section, we explored 
one way to normalize R s , namely by analyzing the ropelength L c of the 
inscribed knot Kp. In doing so, we overestimate the minimum rope- 
length. In this section, we present a different normalization that will 
underestimate the ropelength. In the next section, we combine these 
two notions to obtain a reasonable approximation of the minimum ro- 
pelength with relatively few edges. 

The simplest way to determine the ropelength is to sum the length 
of all edges of the knot and divide by R S (P): 
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Figure 8. The raw polygonal L p and inscribed L c rope- 
length for SONO tightened trefoil knot versus the edge 
length with n=2l, 30, 42, 60, 84, 120, 168, 240, 336, 480, 
672, 960, 1342 and 1920. 



L p = [J2 s ij l R ^ P ) =ndl/R s (P) . 

In what follows we shall refer to L p as the raw polygonal ropelength. 
In calculating the ropelength in such a simplistic manner, we implicitly 
assume that the polygonal knot makes the axis of a rope of radius 1, 
which is obviously wrong. Certainly, if the cylindrical segments of the 
rope are not to overlap, the radius should be smaller. This becomes 
clear for large dl. However erratic the raw polygonal length may seem, 
it will prove to be very useful. 

It is interesting to see how the values of the raw polygonal and in- 
scribed ropelength behave in practice as n increases. For the sake of 
convenience, Figure |H1 shows both values plotted versus the edge length 
dl. 

First notice that as the edge length diminishes, the raw polygonal 
and the inscribed ropelengths apparently converge to a common value. 
We call this value the true ropelength, denoted L^. In spite of its 
erratic nature, the raw polygonal ropelength converges faster than the 
inscribed ropelength. The raw polygonal ropelength L p underestimates 
while the inscribed ropelength L c overestimates the ropelength. Can 
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Figure 9. A tight conformation of the (2,59) torus knot 
found by SONO. A typical situation where one piece of 
the rope is winding around another piece is shown in the 
frame. In simulations, the rope is corrugated as shown 
in the inset. 

one find a weighted average that will accurately estimate the ropelength 
with relatively few edges? Below, we present a heuristic reasoning 
which provides the answer. 

4.1. The problem of finding the right ropelength, an analytic 
approach. By looking at SONO tightened knots, see for instance Fig- 
ure El one may notice that we often deal with pieces of a rope winding 
around other pieces of the rope. At smaller n, the rope is strongly 
corrugated. Thus, to be precise, we deal with pieces of a corrugated 
rope winding around other pieces of the corrugated rope. In a model 
situation shown also in the picture, the central piece is straight and 
the other piece winds tightly around it. By analyzing the raw polygo- 
nal and inscribed ropelength functions in this simplified model, we will 
shed some light on their functioning in the description of tight knots. 

To analyze the problem, we consider a still simpler case of a corru- 
gated rope wound toroidally around a corrugated cylinder. See Figure 

uni 

In the latter case, we may perform a rigorous analysis of the ways 
in which the raw polygonal and the inscribed ropelength estimate the 
true value of the ropelength. The latter is known since it is simply the 
length of a smooth torus winding tightly a smooth cylinder. Assuming 
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Figure 10. (a) A piece of a corrugated rope winds 
toroidally around the corrugated cylinder (a straight 
piece of corrugated rope), (b) In the n — > oo limit, we 
deal with a torus winding around a cylinder. 

that both the cylinder and the torus are made of a tube whose radius 
R — 1, the value is Att. 

Now, let us consider the corrugated torus on the corrugated cylinder. 
Both of them can be seen as unions of spheres. We assume that all the 
spheres have radius R — 1. The distance between consecutive spheres 
is dl. Obviously, for an arbitrary dl there are problems with a clean 
closing of the torus. In what follows we shall concentrate on a short 
piece of the torus, so this is not problematic. 

The best way to wind a piece of a corrugated rope around a straight 
piece of the same rope is to wind it inside the groove between two 
consecutive spheres. Figure [TUl illustrates the situation. Our aim is 
to find the raw polygonal length and the inscribed ropelength of a 
single piece of the rope and compare them with the true length of 
the analogous piece of the torus wound tightly around the cylinder. 
The relations we find here should apply to the situation we face in 
calculating the ropelength of the SONO tightened knots. 

Let us look at Figure [TT] Assume the distance between the spheres, 
i.e. the polygon edge length, is dl. Since the corrugated rope is wound 
within a groove, its spheres are at a distance r < 2 from the center O 
of the axis of the corrugated cylinder. We get 

r = y/A - dl 2 /4 . 

The distance is different when the corrugated rope is wound not 
within the groove, but on the hill of the corrugated cylinder. Here 
r = 2. In real situations, such as this presented in Figure the 




20 



JUSTYNA BARANSKA, PIOTR PIERANSKI AND ERIC J. RAWDON 




Figure 1 1 . The geometry of the corrugated rope wound 
within a groove of a straight piece of the same rope. The 
value of R c is the radius of the circles at which consecu- 
tive spheres intersect. Thus, it is also the radius of the 
smooth rope which can be safely placed inside the cor- 
rugated rope. Its axis runs along the inscribed arcs. 




Figure 12. The bead of the winding rope can be lo- 
cated anywhere on the circular hill between two consec- 
utive grooves. 



corrugated rope runs both in the grooves and on the hills. It is thus 
reasonable to assume an average r. It would be an error to calculate 
the average as an arithmetic average, since the shape of the hills is not 
saw-tooth. The hills have the form of circular arcs. 

The position of the bead can be parametrized by the angle ip, whose 
maximum (p ma x value is given by 

Pmax = arcsin 

See Figure The average value of r is found by integration: 

r av = / 2cos{<f)d(f = - . 

<f 'max Jo 2arcsm(j) 
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Let us return now to Figure [TTJ The angle a at which the edge is 
seen from the center O is given by 

'dV 



a 



2 arcsin (arcsin ( 



4 



The middle point of the edge is found at the distance h from O, 
where 

dl Jarcsm(j) 2 - 1 

h = — . 

2 

The length of the inscribed arc joining the middle points of two 
consecutive edges is 



dL, I ( dl x 



dX = dl arcsin (arcsin ( — )) \ arcsin , 

Since the arc makes the axis of the smooth rope of radius 

Rc=y/1- dP/A , 
the normalized length dL c = dA/R c of the rope segment is 



2 dl arcsin (arcsin (*)) J arcsin( *) 2 — 1 

dL c = , \ ■ 

V4 - dl 2 

The raw polygonal length of the segment is simply 

dL p = dl . 

Now, let us ask the basic question: what is the true ropelength dL 
of the segment? It is the length of a piece of a smooth torus (wound 
tightly on the cylinder of unit radius) seen at angle a from the center 
O: 

dL = 4 arcsin ( arcsin ( — 



(arcsin 



The relationship between the true length dL, its raw polygonal length 
dL p , and its inscribed length dL c estimations should be similar to those 
observed for the true L^, raw polygonal L p , and inscribed L c rope- 
lengths found for knots tightened by SONO. Thus, let us see how 
the raw polygonal dL p and inscribed rope dL c approximations differ 
from the true value dL. To get a clear quantitative estimation of er- 
rors that we make using dL p and dL c , we plot the relative deviations 
(dLp — dL) j dL and (dL c — dL) / dL of the raw polygonal and inscribed 
segment length from its true dL value. See Figure ITH1 
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Figure 13. The relative deviations of the raw polygonal 
dL p and inscribed dL c segment lengths from the true 
value dL. The value dL a is the relative deviation of the 
weighted average dL a = (4dL p + dL c )/5. Compare the 
plots with the plots shown in Figure ED 




Figure 14. Weight a versus dl. 



By looking at the picture, one can clearly see that dL p underesti- 
mates while dL c overestimates the length of the segment. Perhaps an 
appropriately weighted average of the values could provide a better 
estimation of the true length. To check this, we solve the equation 

(adLp + dL c ) /(a + 1) = dL . 

Its solution is 

( dZ^arcsinCf )~ 2 -l\ . , . , dl ss 

2 2 ; arcsm arcsm ^ 

V V^dP J v u " 

a = — j, . 

dl — 4 arcsin(arcsin( j)) 

The functional dependence of the weight a on the edge length dl looks 
rather complex, but its plot versus dl reveals that the dependence is 
very weak - at dl = 1, its value is close to 5, but as dl diminishes, it 
converges quickly to 4. See Figure HU 
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Table 1. The raw polygonal, inscribed and weighted 
average ropelength of the trefoil knot tightened by 
SONO. 



n 


dl 


Lp 


L c 


L a = (AL P + L c )/5 


120 


0.27243068 


32.69168 


32.94095 


32.7415 


240 


0.13637831 


32.73079 


32.79243 


32.7431 


480 


0.06820844 


32.74005 


32.75542 


32.7431 


960 


0.03410665 


32.74238 


32.74622 


32.7431 


1920 


0.01705362 


32.74295 


32.74391 


32.7431 



Table 2. The raw polygonal, inscribed and weighted 
average ropelength of the 5i knot tightened by SONO. 



n 


dl 


Lp 


L c 


L a = {4L P + L c )/5 


180 


0.26182883 


47.12919 


47.45869 


47.20 


360 


0.13105733 


47.18064 


47.26225 


47.20 


720 


0.06554618 


47.19325 


47.21360 


47.20 


1440 


0.03277605 


47.19751 


47.20261 


47.20 



By looking at the plot, one may suspect that the expansion of a(dl) 
should have zero order term equal to 4 and much smaller higher order 
terms. We carried out the expansion to get: 

Since the relations between L p , L c and L^ in tight knots should 
be similar to the relations between dL p , dL c and dL, we arrive to the 
conclusion that the weighted average 

L a = (AL P + L c )/5 

should be almost independent of the edge length dl and, thus, it should 
provide a good estimate of the L^ value. To check the hypothesis, 
we perform a series of tests. First of all, we calculate the weighted 
average for the trefoil knot, whose L p and L c are plotted in Figure 
Hoi Looking at the plot of L a , one clearly sees that, as expected, it is 
almost independent of dl. 

Similar tests have been performed for the next torus knots: 5i, 7i 
and 9i. 

The applicability of the weighted average ropelength to other knots 
requires further analysis. 
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32,70 



0,30 



FIGURE 15. Raw polygonal L p , inscribed L c , and 
weighted average ropelength L a for SONO tightened tre- 
foil knot versus the edge length for n— 240, 480, 960 
and 1920. The data were fitted with with second order 
polynomial curves. 

Table 3. The raw polygonal, inscribed and weighted 
average ropelength of the 7i knot tightened by SONO. 



n 


dl 


L p 


L c 


L a = (4L P + L c )/5 


154 


0.39748 


61.21150 


62.21861 


61.41 


308 


0.19921 


61.35573 


61.60250 


61.40 


616 


0.09966 


61.39211 


61.45352 


61.40 


1232 


0.04984 


61.40131 


61.41663 


61.40 



TABLE 4. The raw polygonal, inscribed and weighted 
average ropelength of the 9i knot tightened by SONO. 



n 


dl 


L p 


L c 


L a = (4L P + L c )/5 


188 


0.40140 


75.46436 


76.72711 


75.72 


376 


0.20118 


75.64443 


75.95430 


75.71 


752 


0.10066 


75.69522 


75.77243 


75.71 


1504 


0.05034 


75.70902 


75.72832 


75.71 
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5. Discussion 

Numerical experiments carried out by the SONO algorithm provide 
us with polygonal knots which can be seen as skeletons of tight knots 
tied on a corrugated rope. The vertices of the polygonal knots are 
the centers of the spherical segments of the corrugated tube. As we 
demonstrated, one can place a smooth tube of a smaller radius inside 
the corrugated tube to obtain a smooth knot tied within a smooth 
tube. We know that L C) the bound on its ropelength, is larger than 
the ropelength of the ideal knot. In observing the ropelength of the 
inscribed knots, we see that it converges with the increasing number of 
vertices to a value L^, which can be seen as an upper bound for the 
ropelength of an ideal knot. Determining the value of Loo is a subtle 
problem and shall not be discussed here; however, by combining the 
results we obtained for finite n with the inscribing algorithm, we are 
able to find an upper bound for the ropelength of a few ideal knots. 
For the trefoil knot, the bound equals 32.744, for the 5i knot it equals 
47.203, for the 7i knot it equals 61.417 and for the 9i knot it is 75.728. 
These are the most precise estimations of the upper bounds obtained 
so far. 

As said above, the L c values we find analyzing the polygonal knots 
delivered by SONO can be seen as upper bounds for the ropelength 
of their ideal smooth conformations. But what about the weighted 
average values L a ? To see their utility we must adopt the less rigor- 
ous point of view of an experimental physicist. The coordinates of the 
polygonal knots delivered by SONO can be seen as results of measure- 
ments carried out on knots tied on the corrugated rope. Analyzing the 
experimental data, one finds that the knot length calculated according 
to the raw polygonal length formula systematically increases with n 
while the length calculated according to the inscribed arcs formula sys- 
tematically decreases with n. Plotting the values together, one notices 
that they converge to a common value: L c from below, L p from above. 
One also finds that, located between the two plots, the plot of their 
weighted average is almost flat and horizontal. By looking at the L a 
values presented in the tables, one can clearly see that L a provides a 
good estimate of the ropelength, even at a small number of vertices. 
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